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Abstract 

'^ ' We investigate the addition of stiffness to the lattice model of hydrogen-bonded polymers 

Q ' in two and three dimensions. We find that, in contrast to polymers that interact via a 

I ^1 ' homogeneous short-range interaction, the collapse transition is unchanged by any amount of 
stiffness: this supports the physical argument that hydrogen bonding already introduces an 

^ ' effective stiffness. Contrary to possible physical arguments, favouring bends in the polymer 

rz ' docs not return the model's behaviour to that comparable to the semi-flexible homogeneous 

OO ' interaction model, where the canonical 0-point occurs for a range of parameter values. In 



fact, for sufficiently large bending energies the crystal phase disappears altogether, and no 
phase transition of any type occurs. We also compare the order-disorder transition from the 
00 , globule phase to crystalline phase in the semi-flexible homogeneous interaction model to that 

for the fully-flexible hybrid model with both hydrogen and non-hydrogen like interactions. 



We show that these phase transitions are of the same type and are a novel polymer critical 
phenomena in two dimensions. That is, it is confirmed that in two dimensions this transition 
Cu I is second-order, unlike in three dimensions where it is known to be first order. We also 

estimate the crossover exponent and show that there is a divergent specific heat, finding 
(p — 0.7(1) or equivalcntly a = 0.6(2). This is therefore different from the 9 transition, for 
which a = —1/3. 
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1 Introduction 

When modelling a single polymer in dilute solution by a lattice self-avoiding walk [H [H E] the 
effects of solvent mediated intra-polymer interaction are often included by assigning an energy 
to each non-consecutive pair of monomers lying on the neighbouring lattice sites: the interaction 
is homogeneous in that in doesn't depend on shape of the two parts of the polymer containing 
the interacting monomers. This is the well-studied Interacting self-avoiding walk (IS AW) model 
which is the standard model of polymer collapse using self-avoiding walks. If the energy is 
repulsive the polymer behaves as a swollen chain (the so-called excluded- volume state) regardless 
of temperature and one says that it is in a good solvent. When the energy is attractive, and the 
temperature is low enough, the chain becomes a rather more compact globule [2lll]j reminiscent 
of a liquid droplet: this is also known as the poor solvent situation. The transition point between 
these two phases is called the 0-point; it is a well studied continuous phase transition (see [5] and 
references therein). However, the mapping of monomers onto lattice sites ignores the natural 
rigidity of real polymers. An energy for bends in the self-avoiding walk can be introduced to take 
account of this feature. The model of semi-flexible polymers has been investigated mostly in 
three dimensions [6l[71[8]. In particular, Bastolla and Grassberger [6] investigated semi-flexible 
interacting self-avoiding walks (semiflexible IS AW) on the cubic lattice, which interact via all 
nearest-neighbours, as in the 0-point model, and included the bending energy. They showed that 
when there is a strong energetic preference for straight segments, this model undergoes a single 
first-order transition from the excluded-volume high-temperature state to a crystalline state. 
Intriguingly, if there is only a weak preference for straight segments, the polymer undergoes two 
phase transitions: on lowering the temperature the polymers undergoes the ^-point transition 
to the liquid globule followed at a lower temperature by a first-order transition to the frozen 
crystalline phase. In two dimensions the transition between the globule and the frozen state has 
been studied in Hamiltonian walks, and there it seems to be a continuous one [H]. 

The modelling of polymers in solution changes as soon as we want to describe any biological 
system (e.g. proteins), in which the hydrogen bonding plays an important role [10]. One of the 
main features of the bonding is that the interacting residua lie on a partially straight segments 
of the chain. Hydrogen-like bonding was first modelled on the cubic and square lattices using 
Hamiltonian paths by Bascle et al. jll| . A monomer acquires a hydrogen-like bond with its 
(non-consecutive) nearest neighbour if both of them lie on straight sections of the chain. Note 
that the identification of a single contact of this type with a single hydrogen bond is only valid 
if fully-flexible polymers are considered; otherwise the contact represents an agglomeration of 
such bonds. The interacting self-avoiding walk modified to have only such interactions will be 
referred to as the hydrogen-like bonding model, or rather IHB model. The IHB model was 
studied in mean-field approximation [11] and a first-order transition from a high-temperature 



excluded-volunie (swollen) phase to a quasi-frozen solid-like phase was found in both two and 
three dimensions. Hence this would indicate that it is a different type of transition from the 
0-point. Note also that the low temperature phase was found to be anisotropic whereas the 
collapsed globule of the standard 6'-point model is isotropic. The IHB model on the square 
lattice was studied directly by Foster and Seno by means of the transfer matrix method [12] and 
by Krawczyk et al. [13] on both the square and cubic lattice using a Monte Carlo method. In 
both of these studies a first-order transition was found between an excluded-volume (swollen- 
coil) state and an anisotropic ordered compact phase in two and in three dimensions, again in 
opposition to the 0-point [2] . 

However, the IHB model was recently extended [TJ] to a hybrid model {IHB-INH) that 
includes both the hydrogen-like bond interactions and non-hydrogen like bond interactions, 
with separate energy parameters. When the non-hydrogen bonding energy is set to zero the 
IHB model is recovered. If both energies are set to be the same then the ISAW without stiffness 
is recovered. For large values of the ratio of the interaction strength of hydrogen-bonds to non- 
hydrogen bonds, a polymer will undergo a single first-order phase transition from a swollen coil 
at high temperatures to a folded crystalline state at low temperatures. On the other hand, for 
any ratio of these interaction energies less than or equal to one there is a single 0-like transition 
from a swollen coil to a liquid droplet-like globular phase. Importantly, for intermediate ratios 
two transitions can occur, so that the polymer first undergoes a ^-like transition on lowering the 
temperature, followed by a second transition to the crystalline state. In three dimensions it was 
found that this second transition is first order, while in two dimensions it is probably second 
order with a divergent specific heat. In other words, at least in three dimensions, by adding an 
energy to both the hydrogen-like and non-hydrogen like interactions a phase diagram similar to 
the one for the semi-flexible ISAW is obtained. 

Various issues then arise. Firstly it is worth studying the addition of stiffness to the IHB 
model since stiffness clearly affects the type of phases that occur in the ISAW model. Moreover, 
one could argue that since in the IHB model interactions occur only between straight segments 
of the walk, an effective stiffness has already been introduced: one could then go further and 
argue that favouring bends may result in behaviour like that in the ISAW model. Secondly, it 
is worth investigating the semi- flexible ISAW model in two dimensions to check if the similar- 
ity of the IHB-INH phase diagram to the semi-flexible ISAW extends to that dimension. In 
particular, given the existence of the same three phases (as we shall find) whether the globule- 
crystal transition is of the same second order type as in the IHB-INH model. Finally the more 
general relationship between these three potentials (nearest-neighbour interaction, hydrogen-like 
bonding and stiffness) on the phase structure of the model polymers is worth pursuing. Hence, 
in this paper we investigate via Monte Carlo simulations the effect of adding stiffness to the 
interacting hydrogen bonding model and compare this to adding stiffness to the canonical ISAW 



model. Intriguingly we find that favouring stiffness does not change the single transition found 
in the IHB model without stiffness. However, for a sufficiently large bending energy favouring 
bends the transition temperature is found to go to zero, and for larger ratios of bending energy 
to hydrogen bonding energy no phase transition occurs. We then compare the globule-crystal 
phase transition in the semi- flexible ISA W model on the square lattice (on the cubic lattice they 
are both first-order) to that of the IHB-INH model and show that the exponents are most likely 
the same: we also estimate these exponents. 

In general we have considered various restrictions of a model of three parameters where 
stiffness is added to the hydrogen-bond interactions and non-hydrogen-bond interactions. In our 
conclusions we argue that only three phases exist in the larger parameter space: an excluded- 
volume dominated state (where the polymer is 'swollen'), a disordered globular state where the 
polymer is in a condensed liquid-like drop and a crystalline state. 

The paper is organised as follows. In Section 2 we explain more carefully details of the 
models considered. In Section 3 we consider the semi-flexible interacting hydrogen bond model. 
In Section 4 we consider the semi-flexible version of the canonical interacting SAW model and 
we compare our results to those from a model that has different energies for hydrogen and non- 
hydrogen nearest-neighbour interactions. In Section 5 we investigate a semi-flexible interacting 
polymer where there are no hydrogen bonds. We end with some conclusions about the most 
general model deflned in Section 2. 

2 Definitions 

Let us deflne a general model that contains each of the models considered as sub-cases via 
restricting parameters. 

We begin with a self-avoiding walk on the square and simple cubic lattices. The walk consists 
of a sequence of occupied lattice sites joined by steps of the walk. A walk of n steps occupies 
n + 1 sites. Consider the sites of the lattice occupied by the walk. When two sites of the walk 
are adjacent on the lattice and not consecutive along the walk, so as not to be joined by a step 
of the walk, we refer to this pair of sites as a nearest-neighbour contact. Additionally, let us 
refer to two consecutive steps that follow the same lattice direction as a stiff step-pair, so that 
there are three consecutive occupied sites along a line on the lattice, and the site in the centre 
of this trio to be a stiffness site (see Figure d]) . 

Now partition the possible types of contact into two classes: when they occur between 
stiffness sites then we refer to these as hydrogen-bond contacts, and all others are non-hydrogen- 
bond contacts. In Figure[2]the partition of the types of nearest-neighbour contacts for the square 
lattice is shown. In Figure [3] the hydrogen-bond contacts (only) are shown for the cubic lattice. 

We now add an energy to the self-avoiding which consists of three contributions: an energy 
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Figure 1: Stiffness sites and non stiffness sites on tlie square lattice 
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Figure 2: The partition of the types of nearest-neighbour contacts into hydrogen-bonds and 
non-hydrogen-bonds for the square lattice is shown. Rotations of these are also possible. 
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Figure 3: Hydrogen-bond contacts for the simple cubic lattices. 

for each hydrogen-bond contact —£hh^ a-n energy of non- hydrogen-bond contact —Enh and an 
energy for stiff step-pairs —ess. The total energy of a walk configuration (^„ of n steps is 

Eni^P-n) = -mhi,{ipn) Ehb - rrinhifn) Enh " rUssifn) Egs, (2.1) 

where ruhb denotes the number of of hydrogen-bond contacts, runh denotes the the number of 



non-hydrogen-bond contacts, and rriss denotes the number of stiffness sites in the walk configu- 
ration (see Figured]). 
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Figure 4: A sample walk configuration with the various type of interaction highlighted for the 
square lattice. 

We will denote the total number of all nearest-neighbours interactions as rriis, where it is 
equal to the sum of the number of the two types of interaction considered in our full model, that 
is rriis = runh + mhb- 

The inverse temperature is denoted as /3 = l/ksT, where ks is the Boltzmann constant and 
T the absolute temperature. We define for convenience Phh = /3e/i6, Pnh = f^^nh and (3ss = P^ss- 
The partition function is then given by 



ZniPhb, Pnh, Ps 
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(2.2) 
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with Cn,mhb,m„h,mss ^^^ density of states. Canonical averages are calculated with respect to this 
density of states. 

Our results are for the following models: 

• the semi-flexible interacting hydrogen-bonding model (semi-flexible IHB model) where 

Pnh = 0, 

• the semi-flexible interacting non-hydrogen-bonding model (semi-flexible INH model) where 
Phb = 0, 

• and the semi-flexible interacting self-avoiding walk model (semi- flexible IS AW model) 
where Phb = Pnh- 



We remind the reader that the semi- flexible IS AW model on the simple cubic lattice has previ- 
ously been studied by Bastolla and Grassberger [B]. 

We compare our results to those concerning the previously studied [ll] model where different 
energies are assigned to hydrogen-bond contacts and non-hydrogen-bond contacts, though not 
to stiffness sites, namely 

• the Interacting hydrogen-bonding - Interacting non-hydrogen-bonding model (IHB-INH 
model) where Pss = 0. 

All the simulations in this paper use a Monte Carlo technique, known as FlatPERM |15j . 
which is well suited to the study of self-avoiding walks on the simple cubic and square lattices 
with interactions. This technique allows for the estimation of quantities at all values of an 
interaction parameter by the estimation of the appropriate 'density of states': e.g. from equation 
()2.2p the estimation of Cn,mub,m„h,'m,ss would allow the calculation of canonical averages for any 
value of the parameters j3hi, , [3nh and (5ss ■ On the other hand each new parameter increases the 
computational cost by at least a factor of n (being the range of the variable — e.g. rrihb when 
including (3hi,). Therefore so as to obtain data for reasonable length walks it is necessary to 
restrict simulations to one or two interaction parameters. 

3 Semi-flexible IHB model {(5nh = 0) 

We have simulated the semi-flexible IHB model on the square and simple cubic lattice, estimating 
an appropriate density of states Cn,mhh,mss-> so that averages can be performed for all values of 
the parameters [ihb and [igs- We have estimated (m/jh) and (mss), which are directly related to 
the internal energy, and the variances in these averages which are related to the specific heat of 
the model. These 'two parameter' simulations were completed up to n = 128 steps. 

Let us start by considering the square lattice data. As has been previously found [121 [T3] at 
Pss = 0, we find that for any fixed value of (3ss there is strongly growing peak in the variance of 
ruhb at a single value of Phb- This is indicative of a phase transition at some position I5^^{j3ss) 
and as the normalised peak of the variance is growing close to linearly with n irrespective of (5ss 
it is indicative of a first order phase transition. To confirm this we considered the distribution 
of miib at the finite-size transition point /?^^ {j3ss ]n). In Figure [S] we show this distribution at 
/3[J(0.5; 128), /^^^^^ (0.5; 64) and (3ll;^{-0.5; 128), /^^^^VO-S; 64). The deepening bimodal distribu- 
tions reinforce our conclusion that the transition is first order regardless of f3ss ■ 

By considering the scaling of the end-to-end distance and the anisotropy parameter (see 
[13]) we have verified that for (3fib << f^hbif^ssjn) the extended phase exists while for f3hb » 
Phbi(^ss,n) the anisotropic crystal phase is observed. 

In Figure [6] a plot of /3)j^ (/Ss^; 128) is given in the two-dimensional space of Phb, and Pss- 
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Figure 5: Distribution of internal energy for the square lattice semi-flexible IHB model at the 
estimated transition point j3f^j^{f5ss;n) when Pss = 0.5 (left), and Pss = —0.5 (right), using 
n = 128 and n = 64. 

Next to that in Figure [6] the same curve is plotted in the three-dimensional space of Phb, Pnh 
and Pss- 

The results for the cubic lattice are completely analogous and the finite size phase diagram 
is given in Figure [71 

On both lattices we therefore have found that the addition of positive stiffness to the inter- 
acting hydrogen-like bond model leaves the single phase transition from a swollen phase at high 
temperatures to a crystalline phase at low temperatures unchanged. When favouring bends, 
so that Ess and Pss ai'e negative, then when any transition occurs it is again of a similar type 
as the fully-flexible model. However there is a range of A = Sss/^hb such that no transition 
occurs. This is relative easy to understand: hydrogen bonds only occur between two stiffness 
sites which are suppressed by large negative values oi Sss- Let us focus on the square lattice. On 
the square lattice for A = Sss/^hb < — 1 no phase transition occurs on lowering the temperature. 
At zero temperatures the ground state for A < —1 is a walk consisting only of bends with energy 
zero (this state should have a positive entropy). For — 1 < A the ground state is one with long 
folds (/3-like sheets) and a bulk energy —n{ehb + ^ss) which is negative so long as Shb > — ^ss- 
This state has zero entropy. The difference in entropy accounts for the apparent shift of the 
asymptote of the phase boundary from Phb = —Pss to Phb = —Pss + c in Figure [3 



4 Semi-Flexible IS AW (/?„/, = f^hb) 

For the semi- flexible IS AW model we focus our attention on the square lattice which has not 
previously been investigated. On the square lattice we performed simulations for n = 128 for two 
parameters Pis and Pss as well as for the one parameter Pss for constant Pis = 0.7 for lengths up 
to n = 512. The line Pis = 0.7 was chosen so as to focus on the transition between two collapsed 
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Figure 6: Plot of the finite size phase boundary for the semi-flexible IHB model on the square 
lattice (left). The plane of the parameters of the semi- flexible IHB model in the more general 
three parameter space: the two phases are denoted by a dark (red) shading for the crystalline 
phase and a mid-density (green) shading for the extended swollen polymer phase, (right) 



phases: the collapsed-globule and the crystal. 

Let us begin with the two parameter simulations. By considering the maximum eigenvalue 
of the matrix of fluctuations of Tiiis and Trigs we have mapped out a finite sized phase boundary, 
see Figure El 

The boundaries clearly divide the phase space into three phases. Once again, by considering 
the scaling of various quantities such as the end-to-end distance at fixed points deep within 
each suspected phase we are satisfied that the three phases are the same as in the IHB-INH 
model dU on the square lattice: a swollen phase with u = 3/4, and two collapsed phases where 
u = 1/2. In one phase the typical configurations are clearly anisotropic, looking like folded 
/3-sheets, indicating that it is a crystalline phase. Hence the phase structure and phase diagram 
is similar to the three-dimensional case. We have attempted to locate the triple point which 
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Figure 7: Plot of the finite size phase boundary for the semi-flexible IHB model on the cubic 
lattice 

seems to around {(3ss,Pis) = (2.0,0.3). However using different methods has resulted in quite 
different estimates so we do not propose any error estimate on these values. 

For convenience and comparison the corresponding finite size phase boundary diagrams for 
the IHB-INH model studied by Krawczyk et al. |14j are given in Figure [9] . 

Since the extended-globule 9 transition has been well studied we have focussed on the two 
other transitions: extended-crystal and globule-crystal. Both these transitions are first order on 
the cubic lattice [6]. In the IHB-INH model [H] the extended-crystal was also first order on 
both the square and cubic lattice. So firstly let us consider this transition in the semi-flexible 
IS AW model on the square lattice. By considering the scaling of the fluctuations in tti^s and 
the distribution of russ we confirm that the first order nature of this transition. Figure [TOl 
demonstrates that the data is consistent with first order scaling at the transition point of the 
swollen to crystal phases. 

In Figure [TT] the maximum of fluctuations for the transition between the globule and crys- 
talline phases for the semi-flexible and IHB-INH models are plotted together. The semi-flexible 
data is taken from the one parameter runs with Pis = 0.7 at lengths up to n = 512 and the 
IHB-INH data from simulations at Pnh = 1-0 at lengths up to 256. While the data still has 
some corrections to scaling present the divergence of the peak of the fluctuations seem to be 
controlled by the same value of exponent. Ignoring corrections to scaling on the length range 128 
to 256 would give us an estimate of the exponent controlling this divergence for the IHB-INH 
and semi-flexible ISAW models respectively, of 0.41(2) and 0.43(2). 

However, by being conservative with taking into account the corrections to scaling for the 
largest lengths we have estimated the exponent to be 0.4(2). Now, this exponent is expected to 
be a(f) where a is the specific heat exponent of the thermodynamic limit transition and (j) is the 
crossover exponent. Additionally it is usually expected that 2 — a = 1/0 so we have a = 0.6(2) 
and (j) = 0.7(1). 
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Figure 8: Plot of the finite size phase boundary for the semi- flexible IS AW model on the square 
lattice (top). The plane of the parameters of the semi- flexible IS AW model in the more general 
three parameter space: the three phases are denoted by a light (blue) shading for the globular 
phase, a mid-density (green) shading for the extended swollen polymer phase and a dark (red) 
shading for the crystalline phase, (bottom) 

5 Semi-flexible INH model {/3hb = 0) 

Given that we have investigated the effect of stiffness on the IHB model and previously [14j 
considered the IHB-INH model [Hj it was desirable to consider the effect of stiffness on the 
INH model for completeness. 

When (3ss = the INH model behaves exactly as the IS AW with a single collapse transition 
from the extended phase at high temperatures to the globular collapsed phase at low tempera- 
tures. While it is more difficult to detect the 9 point in two dimensions since the specific heat 
exponent q < it signature can still be seen in the specific heat data. We have determined a 
finite size phase boundary (3!^j^{(3ss] n = 128) that can be seen in Figure [T2l We immediately see 
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Figure 9: Plot of the finite size phase boundary for the fully- flexible IHB-INH model on the 
square lattice (top). The plane of the parameters of the semi- flexible INH model in the more 
general three parameter space: the three phases are denoted by a light (blue) shading for the 
globular phase, a mid-density (green) shading for the extended swollen polymer phase and a 
dark (red) shading for the crystalline phase (bottom). 

that while the INH model behaves similarly to the IS AW model the addition of stiffness does 
not produce the same effect. At most one phase transition is found for any positive or negative 
stiffness energy. 

On both lattices we therefore have found that the addition of negative stiffness (effectively en- 
couraging bends) to the interacting non-hydrogen-bond model leaves the single phase transition 
from a swollen phase at high temperatures to a globular phase at low temperatures unchanged. 
This can be simply understood by noting that most of the non-hydrogen-bond contacts occur 
between bends in the walk anyway. The transition temperature goes to zero as the stiffness 
energy goes to negative infinity. 

When favouring straight segments, so that e^s and [iss are positive, then when any transition 
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Figure 10: Plots of the fluctuations per monomer divided by N in the number of ISAW contacts 
{rriis) at Pss = 2.0 for 2d with the horizontal axis scaled as {{Pis — l3l^)N). We have used 
f5\f = 0.27. Shown are lengths 64, 128 and 256. 
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Figure 11: The maximum of fluctuations for the transition between the globule and crystalline 
phases for two models: semi-flexible ISAW (a) and the IHB-INH (b) models. (The first curve 
(a) is shifted to show them next to each other.) 

occurs it is again of a similar type as the fully- flexible model (that is 0-point-like). However, at 
least for A = Ess/^nh large enough no phase transition occurs on lowering the temperature. 

A similar argument to the one in Section 3 seems to hold: non-hydrogen-bond contacts do 
not occur between stiffness sites which are favoured by large positive values of Ess- We note that 
the zero temperature states are pathological and there exits zero temperature phases transitions 
(rod-coil and rod-globule). Nevertheless it seems that the change in the zero-temperature state 
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on varying the parameters still accords with the position of the finite temperature swollen- 
globule phase boundary. At zero temperatures the ground state for A > 2 is a walk consisting 
only of straight segments with energy — ne^s- For < A < 2 the ground state consists of long 
zig-zag paths with each 'zig' and each 'zag' being made up of two steps (in this way each straight 
segment is always adjacent to two bends) next to each other which have one non- hydrogen bond 
per step and one stiffness parameter per two steps: the energy of this state is —n£hb — ^£ss- 
These two states cross energies at A = 2: this seems to explain the asymptote of the phase 
boundary in Figure [12] for large positive Ess- 
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Figure 12: Plot of the finite size phase boundary for the semi-flexible INH model on the square 
lattice (top). The plane of the parameters of the semi- flexible INH model in the more general 
three parameter space: the two phases are denoted by a light (blue) shading for the globular 
phase and a mid-density(green) shading for the extended swollen polymer phase, (bottom) 
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6 Conclusions 

We have investigated the effect of increasing stiffness and also enhancing bends on the lattice 
model of hydrogen-bonded polymers. We have found that in both cases if there is a phase 
transition it is unchanged from the fully-flexible model: namely a first order phase transitions 
occurs in both two and three dimensions. We also argue that if bending is sufficiently enhanced 
no phase transition occurs at all. This is in contrast to the effect of adding stiffness to the 
canonical model of self-interacting polymers where adding stiffness results in three phases: a 
high temperature excluded-volume dominated "swollen" phase, a liquid-like globule phase and 
an anisotropic solid- like polymer crystal phase. We have investigated this semi- flexible IS AW 
problem in two dimensions and shown that these three phases exist as they had previously been 
shown to exist in three dimensions. We have investigated the globule-crystal transition on the 
square lattice more closely and found that unlike three dimensions where it is first order but 
like another recently studied model, extending the hydrogen bonding model by the addition 
of non-hydrogen bond interactions, the transition is second order with specific heat exponent 
a = 0.6(2). 

Putting together all the information at hand it is likely that in the three-dimensional phase 
space of hydrogen-bond, non-hydrogen-bond nearest neighbour interactions and stiffness only 
the three phases already studied occur. In Figure [13] all the phase boundaries found when the 
energies are all positive are illustrated. One can the infer that for Phb and Pnh small no matter 
what the value of (3ss the extended phase exists. Also, one can infer that for large f3nh the 
globular phase exists and for large (3nh the crystal phase exists. In this way the partial results 
in the literature can now be understood. 
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